Probing two topological surface bands of Sb 2 Te 3 by spin-polarized photoemission spectroscopy 
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Using high resolution spin- and angle-resolved photoemission spectroscopy, we map the electronic structure 
and spin texture of the surface states of the topological insulator Sb2Te3. In combination with density functional 
calculations (DFT), we directly show that Sb2Te3 exhibits a partially occupied, single spin-Dirac cone around 
the Fermi energy Ef, which is topologically protected. DFT obtains a spin polarization of the occupied Dirac 
cone states of 80-90 %, which is in reasonable agreement with the experimental data after careful background 
subtraction. Furthermore, we observe a strongly spin-orbit split surface band at lower energy. This state is found 
at E — Ef — — 0.8 eV at the T-point, disperses upwards, and disappears at about E — Ef = — 0.4 eV into two 
different bulk bands. Along the F — K direction, the band is located within a spin-orbit gap. According to an 
argument given by Pendry and Gurman in 1975, such a gap must contain a surface state, if it is located away 
from the high symmetry points of the Brillouin zone. Thus, the novel spin-split state is protected by symmetry, 
too. 

PACS numbers: 71.20.Nr, 71.10.Pm, 71.70.Ej, 73.20.At 



I. INTRODUCTION 

Topological insulators (TIs) are a new phase of quantum 
matter giving rise to, e.g., a quantum spin Hall phase with- 
out external magnetic field 0H3- Large spin orbit (SO) in- 
teraction and inversion symmetry lead to nontrivial edge or 
surface states which reside in a bulk energy gap and are pro- 
tected by time reversal symmetry. This new phase is classi- 
fied by a Z 2 topological invariant, which distinguishes it from 
an ordinary insulator (4). In three dimensions, the surface 
states form an odd number of massless Dirac cones exhibit- 
ing a chiral relationship between spin and momentum of the 
electrons [5-9]. Bii_ x Sb x was the first 3D TI to be theoreti- 
cally predicted [5, 6] and experimentally discovered |[T0l[TT1l . 
Subsequent calculations showed that the thermoelectric mate- 
rials Bi 2 Se3, Bi 2 Te3 and Sb 2 Te3 should exhibit even simpler 
TI properties with only one Dirac cone around T lfl2l . After- 
wards, Bi 2 Se3 and Bi 2 Te3 were studied by spin- and angle- 
resolved photoemission spectroscopy (spinARPES) confirm- 
ing the TI properties of the surface state fT3llT8ll . Scan- 
ning tunneling spectroscopy (STS) revealed the absence of 
backscattering, i.e. momentum inversion, for the surface state 
of these materials, which confirms the protective spin chi- 
rality of the Dirac cone fT9U2Tll . While the TI nature of 
Bi 2 Se3 and Bi 2 Te3 is well established, the phase change ma- 
terial Sb 2 Te3 l22l was rarely probed. Recent ARPES mea- 
surements suggest that the Fermi level, E F , of single-crystal 
Sb 2 Te3 is within the bulk valence band probably due to hole 
doping fT6ll . In contrast, thin films grown by molecular beam 
epitaxy (MBE) exhibited a Dirac cone around Ep in ARPES 
l23l . However, the spin chirality and the topological nature of 
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the Dirac cone has not been tackled so far. 

Here, we report spin resolved ARPES data in combination 
with densitiy functional theory (DFT) calculations on single 
crystal Sb 2 Te3 showing the spin-polarized nature of the Dirac 
cone. In addition, we found another surface state at lower en- 
ergy exhibiting strong Rashba-type spin- splitting. This split- 
ting described by the Rashba parameter a ~ 1.4 eVA is larger 
than for Au(lll) (HIES or Bi(lll) I26ti28l but lower than 
in Bi-based surface alloys E3 EOj. DFT calculations repro- 
duce this state and reveal that it is located in a narrow gap 
in the T — K direction that originates from SO interaction. 
This observation is in line with the analytical prediction from 
1975 [31], that in the case of a SO gap away from high sym- 
metry points of the Brillouin zone, there must be, at least, one 
surface state present. This state is, thus, the second surface 
state of Sb 2 Te3, which is protected by symmetry. 



II. EXPERIMENT 

Spin- and angle-resolved photoemission (spinARPES) ex- 
periments have been performed at 300 K with electron analyz- 
ers Scienta R8000 and SPECS PHOIBOS 150 using linearly 
polarized synchrotron radiation from the beamlines UE112- 
PGM-1 and UE112-lowE-PGM2 at BESSY II. Energy and 
momentum resolution of the Scienta R8000 analyzer are 
20meV and 0.6 % of the surface Brillouin zone, respectively. 
Spin resolution is achieved with a Rice University Mott po- 
larimeter operated at 26 kV (Sherman function S e ^ = 0.16) 
and capable of recording the two orthogonal spin directions in 
the surface plane of the sample. The energy and momentum 
resolution of the spin resolving apparatus are lOOmeV and 
4 % of the surface Brillouin zone, respectively. Scanning tun- 
neling microscopy (STM) is performed with a modified Omi- 
cron STM at 300 K. For all measurements, the Sb 2 Te3 single 
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FIG. 1. (color online) (a) STM image of cleaved Sb 2 Te 3 (0001); V = 
-IV, I = 0.8 nA, 375 x 375 nm 2 . Inset: Line profile showing the 
step which corresponds to the height of one QL. (b) ARPES data of 
Sb2Te3(0001) along T — K at an incident photon energy hv — 50 eV; 
Dirac cone is marked. 



crystal is cleaved in ultra high vacuum (UHV) at a base pres- 
sure of 1 • 10 -10 mbar. Large terraces with widths of several 
100 nm separated by step edges of 1 nm height result as veri- 
fied by STM (Fig.[TJa)). Figure[TJb) reveals that ARPES mea- 
surements exhibit distinct bands after cleavage, however with 
a relatively strong dependance on photon energy (see below). 

The calculations are performed within the generalized gra- 
dient approximation l32li to DFT, employing the full-potential 
linearized augmented planewave method as implemented in 
the Fleur code l33l . SO coupling is included in a non- 
perturbative manner (34). Based on the optimized bulk lattice 
parameters, the surfaces are simulated by films of a thickness 
of six quintuple layers (QLs) embedded in vacuum. 



III. GENERAL CHARACTERIZATION 

The sketch in Fig. |5Jc) shows the crystal structure of 
Sb2Te3 consisting of consecutive QLs with stacking sequence 
Te(l)-Sb-Te(2)-Sb-Te(l). The different numbers in parenthe- 
ses mark the different environments of the Te layers. The cou- 
pling within a QL is strong, whereas the interaction between 
two QLs is predominantly of van der Waals type JT2) . Conse- 
quently, cleavage leads to a Te terminated (0001) surface with 
hexagonal symmetry as has been verified by low-energy elec- 
tron diffraction. Identically to Bi 2 Se3 and Bi 2 Te3, Sb 2 Te3 has 
inversion symmetry with the layer Te(2) containing the cen- 
ter of inversion. This simplifies the calculation of Z 2 , which 
becomes an analysis of states at the high symmetry points 
only, and, thus, leads to the straightforward identification of 
a strong topological insulator (8lfT2). 

The electronic structure of Sb 2 Te3 is firstly probed by STS 
as shown in Fig. [2ja). STS records the differential tunneling 
conductivity dl/dV which is proportional to the local density 
of states (LDOS) of the sample [ 35 1 . As theoretically pre- 
dicted (T2) . a small band gap of about 200 meV is observed. 
The Fermi level Ep (sample bias V = V) is found close 
to the valence band edge indicating p-type doping (T6J |23) 
in agreement with a recent STM study on MBE grown thin 
films 1361. The band structure from DFT calculations of a six 



QL film including SO interaction is shown in Fig. [2fb). It is 
in good agreement with a recent calculation (37). Small dif- 
ferences to earlier calculations (121 [16) can be traced back to 
the sensitivity of the electronic structure to small changes of 
the geometrical parameters. The spin-polarized states are dis- 
played as colored circles being blue or red for the different 
spin orientations. Only the spin polarization perpendicular to 
the in-plane wave vector of the electrons kn and the surface 
normal is shown. The varying radius marks the absolute value 
of this k-resolved spin density at and above the surface. The 
absolute spin polarization of the states in comparison with 
the values from the experiment is discussed in more detail 
in Fig. [4] Projected bulk bands resulting from a bulk calcu- 
lation are shown as shaded areas Fig. [2jb)- A single Dirac 
cone originating from topological protection is found around 
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FIG. 2. (color online) (a) dl/dV spectrum recorded by STM on 
the cleaved Sb2Te3(0001) surface; V s tab — —0.5 V, / sta b = InA, 
Vmod = 10mV; band edges of the valence band (VB) and conduc- 
tion band (CB) are marked by lines, Fermi level Ef is marked, (b) 
Band structure in V — K direction as calculated by DFT including 
spin-orbit coupling; states resulting from a film calculation are shown 
as circles with the color (blue or red) indicating different spin direc- 
tions and the size of colored circles marking the magnitude of the 
spin density (for absolute spin polarization values see Fig. [4]); shaded 
areas are projected bulk bands originating from a bulk calculation, 
(c) Sketch of the crystal structure of Sb2Te3; one QL is marked with 
different atoms in different colors as indicated, (d) ARPES data of 
the lower Dirac cone along T — K and bulk valence band (BVB) as 
marked; hv — 55 eV; orange, dashed line is a guide to the eye from 
which the Fermi velocity vf is deduced; white dashed lines mark po- 
sitions of the energy distribution curves (EDCs) shown in Fig. [3|a), 
(b) and (c). 
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T with Dirac point at and overlap of the occupied states 
with bulk states. Furthermore, we checked the Z 2 invariant 
with the help of our DFT calculations and found that Sb2Te3 
is indeed topologically non-trivial with an index of (1;000). 
Another bulk band gap exists around E — E F = —400 meV. It 
houses two spin-polarized surface states exhibiting a Rashba- 
type spin splitting AE = a • \kn | with the wave number par- 
allel to the surface kn and the Rashba coefficient a ~ 1.4 
eVA, at least, up to kn ~ 0.05 A -1 . This a is larger than the 
value for Au(l 11) (a = 0.33 eVA) El orBi(lll) (a = 0.55 
eVA) [27], both consisting of heavier atoms, but lower than 
the largest a- values so far found in Bi surface alloys (a — 3.8 
eVA) [29]. 
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IV. DIRAC CONE IN FUNDAMENTAL BAND GAP 
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Figure |2jd) shows the measured energy dispersion around 
the Dirac point along T — K, which is most easily visible at 
hv w 55 eV. The linear dispersion of the Dirac cone with 
a crossing point close to Ey is visible. The latter indicates 
that the position of the surface Fermi level is predominantly 
determined by the Dirac electrons and not by extrinsic dop- 
ing. This deviates from previous ARPES results obtained on 
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FIG. 3. (color online) (a), (b) Spin-resolved energy distribution 
curves (EDCs) for the spin component perpendicular to ky recorded 
at A; || - values as indicated and marked by dashed lines in Fig.[2]cl); dif- 
ferent colors mark different spin directions; hv — 54.5 eV. (c) Spin- 
resolved EDCs for the spin component parallel to ky ; hv — 54.5 eV. 
(d) Resulting spin polarization perpendicular to the two different k\\ 
as marked, (e) Sketch of the lower Dirac cone with spin directions 
marked as deduced from spinARPES and in accordance with DFT. 



FIG. 4. (color online) Expectation values from DFT for the out-of- 
plane and in-plane component of the spin polarization along V — K 
for the Dirac cone and the Rashba- type surface state (Rashba). Upper 
and lower indicates the energies above and below the Dirac point and 
the energetically higher and lower Rashba band, respectively. 



bulk Sb2Te3 (T6), but is in agreement with ARPES data from 
thin films grown by MBE l23ll . It points to a low defect den- 
sity of the investigated crystal [36]. The linear dispersion 
is fitted by E — E F = hvF\kn\ resulting in Fermi veloc- 
ity vf — 3.8 ± 0.2 • 10 5 m/s (dashed, orange line), which 
agrees reasonably with = 3.2 • 10 5 m/s obtained by DFT 
[Fig. [2jb)]. The background of bulk valence bands (BVB) 
found in DFT is also visible in Fig.[2jd). In order to verify the 
chiral spin polarization of the Dirac cone, we use spinARPES 
sensitive to the spin components within the surface plane. 

Figure |3ja), (b) and (c) show the spin-resolved energy dis- 
tribution curves (EDCs) as measured at fey = —0.06 A -1 and 
fey = 0.06 A -1 as marked. The spin component perpendic- 
ular to ky and the surface normal [Fig. [5ta), (b)] exhibits an 
intensity difference between spin up ana spin down compo- 
nent which reverses for the opposite wave number. In con- 
trast, the spin component parallel to ky [S x , Fig.[3jc)] shows 
no spin polarization. This leads to the spin momentum rela- 
tion depicted in Fig. |3je) for the in-plane spin polarization 
which is typical for a topological insulator. The spin is per- 
pendicular to ky and rotates counterclockwise for the lower 
part of the Dirac cone as in the case of Bi 2 Te3, Bi 2 Se3 f\3\- 
021 or the tunable topological insulator BiTKSi-^Se^ l38l . 
The same sense of rotation is found by the DFT calculations 
[Fig. [2jb)]. Figure [3jd) shows the resulting spin polariza- 
tion for the two opposite momenta calculated according to 
P y = (Si - S$)/(Sl + Sj) with fit and S± being the mea- 
sured spin resolved intensities perpendicular to ky . The polar- 
ization of P y ~ 20 % has opposite sign for opposite momenta. 

The spin polarization of the Dirac cone in Sb 2 Te3 is fur- 
ther calculated by DFT and analyzed in terms of an in-plane 
and an out-of-plane component. Only the in-plane compo- 
nent perpendicular to ky is considered, since we did not find 
any spin polarization for the direction parallel to kn. Figure 
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FIG. 5. (color online) (a) ARPES dispersion of the Rashba type surface state (Rashba SS) along T — K direction; hv — 22 eV; band structure of 
the Rashba surface states from DFT is superimposed as blue and red lines; dashed lines mark the position of the EDCs in (b). (b) Spin resolved 
EDCs (points in red and blue for the two spin directions perpendicular to ky ) at different momenta as indicated and marked by dashed lines in 
(a); fits of the peaks are shown as solid lines (red, blue); peak positions are marked by dashed lines and lead to the spin splitting energies AEso 
indicated; (c) calculated spin splitting (Theo) of the Rashba state in comparison with measured spin splitting (Exp), (d) Two-dimensional cut 
through the calculated local density of states for the Rashba state and the lower and upper part of the Dirac cone at k\\ = 0.06 A -1 . 



[4] shows the resulting spin polarization values with respect 
to the wavenumber integrated over the first two atomic lay- 
ers. This area approximately corresponds to the penetration 
depth in the ARPES experiment at a photon energy of 55 eV. 
While P y ~ 1 is found for free Dirac cones by DFT (T8), 
we find a reduced polarization for the lower Dirac cone of 
roughly 80 % near the T point, which increases to about 90 % 
at fey = 0.06 A -1 . This is mostly due to a penetration of 
the Dirac cone states into subsurface layers (Fig. [5jd)) where 
fluctuating electric fields lead to a complex spin texture. In 
contrast, the in-plane polarization of the upper Dirac cone de- 
creases towards higher wavenumbers down to 60%. Note, 
the considerable out-of-plane polarization in V — K direction 
for higher wavenumbers, which is in line with the warping of 
the Dirac cone at higher energies HUES]. This result agrees 
also qualitatively with the calculations of Yazyev et al. l40l 
for Bi 2 Te3 and Bi 2 Se3. However, the polarization values in 
those materials are typically smaller, reflecting the stronger 
spin-orbit entanglement caused by the heavier Bi atom. 

The discrepancy in the in-plane spin polarization between 
calculation (P y ~ 90 %) and experiment (P y ~ 20 %) at fey = 
0.06 A- 1 can be traced back to the finite angular resolution 
of the spin ARPES experiment. Deconvolution from the BVB 
reveals an estimated spin polarization for the lower Dirac cone 
of 80-90 % in experiment (see appendix), which is in a good 
agreement with the DFT result. 



V. RASHBA SPIN-SPLIT SURFACE STATE PROTECTED 
BY SPIN-ORBIT 

Next, we demonstrate that Sb2Te3 exhibits an additional 
spin split surface state originating from SO interactions. Fig- 
ure [5ja) shows ARPES data along T — K recorded at lower 
photon energy, hv = 22 eV, and superimposed by the DFT 
band structure of the Rashba- type band. Based on the ex- 
cellent concurrence, we conclude that the ARPES data at 



hv = 22 eV are dominated by these bands, while the Dirac 
cone is barely visible. We determine the spin splitting of the 
Rashba bands using spin-resolved EDCs for the spin direction 
perpendicular to ky as shown in Fig. [5jb). Therefore, we fit 
the peak in each curve by a Lorentzian function as shown by 
solid lines. This leads to the spin splitting energies AEso 
indicated and plotted for different ky in Fig. J5jc). Reason- 
able agreement between theoretical and experimental AEso 
is found. Moreover, as expected for a Rashba-type spin split- 
ting, the spin direction for upper and lower peak inverts by 
inverting the ky direction. We checked that negligible spin po- 
larization is found parallel to ky which implies that the spin of 
the upper (lower) band rotates clockwise (counter-clockwise) 
with respect to kn in agreement with the DFT results. As 
depicted in Fig. (4j DFT further reveals that the Rashba state 
shows no pronounced out-of-plane polarization within the first 
two atomic layers whereas in the in-plane direction, the dif- 
ferent spin branches are nearly fully spin polarized. From 
fey = 0.05 A -1 to fey = 0.15 A -1 the lower spin branch 
shows a slightly lower in-plane polarization than the upper 
spin branch. This is probably due to the proximity to the bulk 
band. Altough there is a small reduction of the polarization, 
the coupling to the bulk band seems to be rather low. Other- 
wise, the reduction should be more pronounced. Compared to 
the DFT, the experimental spin resolved data reveal a net spin 
polarization of P y ~ 45 % for the upper band and P y ~ 18 
% for the lower band (Fig.[5](b)). Most likely, the overlap of 
the two bands in the experiment as well as its overlap with the 
bulk bands caused by the limited energy and momentum reso- 
lution of the spin detector is responsible for the small num- 
bers. The fact that the peak at higher energy is sharper is 
probably related to its larger separation from the bulk bands 
[as visible in Fig.[2jb)] leading to longer lifetime. 

The charge density of the Rashba-type state, as shown for 
fey = 0.06 A -1 in Fig.pTd) reveals that it has predominantly 
Te p z character and is localized strongly within the Te sur- 
face layer. In contrast, the states of the Dirac cone are more 
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FIG. 6. (color online) (a) ARPES data including the Rashba state at higher |ky |; hv = 54.5 eV; for better visibility, the derivative with respect 
to E is shown; Dirac cone and Rashba-type surface state (Rashba SS) with superimposed E(k\\) curves of the Rashba states (colored lines) 
from DFT calculation are marked. Dashed boxes indicate the same area in (a) and (b) as well as the |k|| | -position of the EDC in (c). The width 
is given by the angular resolution of the experiment, (b) Bulk band structure along T — K calculated with (gray lines) and without (yellow 
lines) SO interaction; black lines mark the band edges with (dashed) and without (solid) SO interaction; with SO interaction, a gap opens near 
the center of the graph around E — — 0.5eV and k\\ — ±0.26 A -1 ; Y^~\ Y^ mark the parity of the two bands at the Y point; (c) Spin 
resolved EDCs measured at the position of the SO gap (|k|| | as marked); hv — 54.5 eV; peak positions as determined from Lorentzian fits 
(solid lines) are indicated by dashed lines; the resulting AEso is marked. 



Sb p z like and are penetrating more strongly into the bulk of 
Sb2Te3. The different penetration might be the reason why the 
Dirac cone is more easily observed at higher photon energy, 
while the Rashba state dominates the spectra at hv — 22 eV. 
The electric field between the surface Te-layer and the subsur- 
face Sb-layer is deduced from the calculated surface core level 
shift to be about 2 • 10 8 V/m having a strong dipolar contribu- 
tion between Te 5_ and Sb 5+ . This surface dipole is probably 
responsible for the relatively large Rashba coefficient, similar 
to the findings in surface alloys 1 29 ] and layered bulk com- 
pounds ATI . 

Differently from the Rashba bands found so far [ 24 , 23, the 
DFT calculations shown in Fig.[2|b) predict that the different 
spin branches disperse into different projected bulk continuum 
bands. Thus, each spin branch connects the upper and the 
lower bulk band surrounding the gap by dispersing from fey = 
-0.28 A" 1 to fey = 0.28 A" 1 . At T, the surface state is spin 
degenerate as requested by time reversal symmetry. 

Figure [6ja) shows the measured band structure close to the 
point where the Rashba bands merge with the bulk bands ac- 
cording to DFT. Indeed, a band moving upwards and a band 
moving downwards are discernible up to about |k||| = 0.27 
A -1 . Fig. |6jc) reveals that a spin splitting of about 81 meV 
is still visible at |k||| = 0.26 A -1 (in comparison to a spin 
splitting of 80meV in the DFT), i.e., close to the point where 
merging of surface bands and bulk bands is obtained in the cal- 
culation. The origin of this remarkable behavior is illustrated 
in Fig.[6jb), where the calculated bulk band structure with and 
without SO interaction is shown. Obviously the SO interac- 
tion opens a gap between the projected bulk states originating 
from a band Y^ near the Fermi level and a lower- lying Y^ 
band, where (+) and (— ) marks the parity of the states at Y. 
The gap is found at fen =0.26 A -1 along the line r — E of the 
bulk band structure (r — K in terms of surface Brillouin zone). 



In such a SO gap, according to a theoretical argument given 
by Pendry and Gurman 1311 , at least one surface state with 
two spin channels must exist. The requirement is that the gap 
is not located at a high symmetry point of the Brillouin zone 
as in our case. Thus, the observed Rashba split surface state, 
which we experimentally resolve within the spin-orbit gap, is 
protected by this gap. Similarly to the topologically protected 
Dirac cone, it connects the lower and the upper bulk bands, 
which are inverted by SO interaction. So far, there has been 
very little experimental proof of such a surface state within a 
SO gap away from a high symmetry point. Feder and Sturm 
l42t report of a spin-orbit generated gap along the symme- 
try line TH in W(001) by means of tight-binding calculation, 
in which a surface state is found experimentally (43) . An- 
other example is a well definded surface state located in the 
r - T direction in Bi(lll) (44l|45]|. First priciple calculations 
by Gonze et al. l46ll resolved a SO gap on this symmetry 
line which is consistent in energy with the measured surface 
state. Differently from these states, our observation reveals a 
spin-splitted surface state with topological character, which, 
moreover exists in parallel to a surface Dirac cone in the fun- 
damental gap. Thus, it adds a distinct example to Pendry and 
Gurman 's criteria. 

Figure [7] shows the band structure of Sb2Te3 from DFT 
without SO interaction. The surface state at the Fermi level 
does not exist in this case, confirming the topological nature 
of this state. The Rashba state becomes a spin degenerate state 
and merges for positive and negative momenta with the same 
bulk band. Thus, it does not connect the upper and the lower 
bulk band, i.e., it loses its protected character. This proves that 
the unconventional behavior towards higher wavenumbers is 
driven by SO coupling in line with Pendry and Gurman 's ar- 
gument l3T1l . 

In summary, spinARPES reveals the spin texture of the 
Dirac cone within the fundamental gap of Sb2Te3, which ro- 
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FIG. 7. (color online) DFT band structure calculated for a slab with- 
out SO coupling. The surface state (red dots) still exists, but is only 
connected to the lower bulk band exhibiting character. 



tates counter-clockwise for the lower part of the Dirac cone. A 
low defect density of the crystal allows to follow this state up 
to the Dirac point. In addition, in accordance with DFT calcu- 
lations, we identified a novel, strongly spin-split Rashba-type 
surface state which is protected by a spin-orbit gap away from 
T and connects an upper and a lower bulk valence band. This 
state is protected by symmetry according to a fundamental cri- 
terion given by Pendry and Gurman in 1975. 
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APPENDIX: EFFECTIVE SPIN POLARIZATION OF THE 
LOWER DIRAC POINT DEDUCED FROM EXPERIMENT 

Whereas calculation predicts an in-plane spin polarization 
of 90 % for the lower part of the Dirac cone at fen = 0.06 A -1 
[Fig.|4|, we only detect a value of P y ~ 20 % in the experi- 
ment. In most experimental works, the reduction of the spin 
polarization is due to extrinsic factors, like the insufficient in- 
strumental resolution in spinARPES measurements especially 
in the case where different states are relatively close to each 
other (TT] [T3. For surface states well separated from the 
bulk states, it has been demonstrated that the spin polariza- 
tion reaches unity |18|. In our case, the unpolarized back- 
ground from the bulk valence band (BVB) states consider- 
ably reduces the spin polarization. Besides this background, 
the spin resolved spectra feature a finite background from the 
Rashba-type surface state visible as increasing intensity to- 
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FIG. 8. (color online) (a), (b) Spin-resolved energy distribution 
curves (EDCs) for the spin component perpendicular to k\\ after sub- 
traction of the Rashba state background (compare with Fig. 2(a) and 
(b)). (c) Corresponding spin polarization as a function of energy. The 
subtraction of the background increases the averaged spin polariza- 
tion from 20 % to 27 %. 



wards higher binding energies in the spin-resolved spectra in 
Fig. 2 (a) and (b). After subtracting this background by using 
a Gaussian function, the spin polarization reaches a value of 
P y ~ 27 % for k\\ = -0.06 A" 1 (see Fig. ^. 

In order to discuss the absolute spin polarization more 
quantitatively, it is necessary to evaluate the finite contribu- 
tion from the spin-degenerate BVB close to the lower Dirac 
cone to the ARPES intensity (47| (Fig. 1(b) and (d)). 

Figure [9] (a) shows a detailed ARPES measurement of the 
Dirac cone recorded with high resolution. The Dirac cone and 
the adjacent BVB are visible. If this data is convoluted with a 
two dimensional Gaussian function having the energy and mo- 
mentum resolution of the spin resolving detector as FWHM 
in the two directions, the Dirac cone part overlaps much more 
strongly with the BVB (Fig. [9](b)). Indeed, a similar broad- 
ening of the structure is found in the ARPES data recorded 
with the spin detector (Fig. [9] (c)), albeit with weaker inten- 
sity due to the lower efficiency of the spin resolving apparatus. 
Constant-energy cuts from the convoluted ARPES data (Fig. 
ge)) and the spinARPES data (Fig. [9] (f)) exhibit a similar 
distribution, too. As a result, we feel enabled to deduce the 
contribution from the BVB out of the high resolution ARPES 
data, where we can clearly distinguish Dirac cone peak and 
BVB background. The quantitative analysis is shown in Fig- 
ure [9] (d) highlighting the contributions of the Dirac cone and 
the BVB as fitting curves using a Lorentzian for the Dirac 
cone states and a tanh for the BVB, respectively. In the case 
of the Dirac cone, the resulting peak is wider than the mo- 
mentum resolution of the high resolution detector. Thus, as- 
suming a Lorentzian function, which takes into account the 
lifetime broadening of the surface state, is justified. The BVB 
contribution is approximated by a tanh function implying a 
BVB of constant density within the Dirac cone as found in the 
DFT calculations and with no bulk bands beyond the Dirac 
cone. The latter implies the additional restriction within the fit 
that the reversal point of the tanh function matches the peak 
position of the Lorentzian. Width and height of the two func- 
tions are used as independent fitting parameters. From the two 
fit curves, shown in Fig. [5Jd), too, we estimate the spectral 
weight of the Dirac cone states with respect to the total spec- 
tral weight within the spinARPES measurement. Therefore, 
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FIG. 9. (color online) (a) High resolution ARPES data of the Dirac 
cone (close-up view from Fig. 1(d)) at an incident photon energy 
hv — 55 eV. (b) data from (a) convoluted with a Gaussian curve 
taking into account the energy resolution (FWHM: 100 meV) and 
momentum resolution (FWHM: 0.09 A -1 ) of the spin detector, (c) 
ARPES data measured with the spin detector, hv — 54.5 eV. (d)— (f) 
Constant-energy cuts through the Dirac cone along the green line in 
the ARPES data aside. Constant energy cuts in (e) and (f) show a 
similar distribution. In (d), the result of the two component fitting to 
the spectrum is shown as marked Dirac cone and BVB (bulk valance 
band). Ratio of the spectral weight of the Dirac cone Wdc to the 
total spectral weight Wtotai is given. 



we determine the area of each fit curve in Fig. [9jd) within the 
width of the angular resolution of the spinARPES experiment 



of 0.09 A -1 around the probed fey-value of 0.06 A -1 . They 
are called VKdc and Wbvb> respectively. We checked that the 
energy resolution of the experiment is not relevant for the an- 
gular broadening. The result is Wdc/(Wdc+W"bvb) = 0.33 
as indicated in Fig. [9fd). The resulting intrinsic spin polariza- 
tion of the Dirac cone in-plane component is then given by 
the experimentally measured spin polarization P y divided by 
Wdc/ Wtotai leading to a value of Pdc — 82 %. The devel- 
opment of the root-mean square error (RMSE) of the fitting, 
as shown e.g. in Fig. [9jd), is then analyzed for fixed differ- 
ent relative heights of the Lorentzian and the tanh function, 
which corresponds to different Pdc- The RMSE is shown as 
a function of Pdc in Fig- 10 revealing that a spin polarization 
of 80-95 % of the Dirac cone states is most likely. Notice the 
nice agreement with the DFT result of Pdc = 90 %. Notice, 
moreover, that the energetic width of the Dirac cone peak in 
spinARPES (Fig. [8]) is similar to the total width of the BVB 
in the DFT calculations shown in Fig. 1(b) of the main text, 
which supports that our estimate of Pdc is reasonable, i.e. the 
BVB contribution leads to a peak of similar width as the Dirac 
cone peak. In turn, high resolution spinARPES experiments 
are required to measure the intrinsic spin polarization of the 
Dirac cone directly, i.e., without relying on any assumption. 




80 90 100 110 

spin polarization (%) 

FIG. 10. (color online) Root-mean- square error (RMSE) for different 
fits of the spectrum of Fig. |9](d) leading to different WDc/W to tai 
and, thus, to different spin polarizations of the Dirac cone states Pdc • 
The RMSE increases by a factor of three for a spin polarization below 
80 %. 
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